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A.  Summary  of  Work  Accomplished 

r~ 

A  new  second- auantiz^t|on  representation  was  derived  f or llow- energy 
positron-hydrogen  scattering?-^ This  work,  which  was  published  in  PhyS .  Rev. 
A  26,  217  (1982)  and  released  as  Technical  Report  QFT-2  of  this  contract, 
served  as  the  foundation  for  our  subsequent  calculation  of  differential 
cross  sections  for  the  positronium  formation  reaction  e+  +  H  -*  Ps  +  H+. 

The  applicability  of  our  representation  for  positron-hydrogen  scattering  to 
the  above-mentioned  reactive  process  stems  from  the  fact  that  our  novel 
second-quantization  representation  treats  all  channels  (both  scattering  and 
reactive)  on  an  equal  footing.  The  calculations  for  the  positronium- 
formation  reaction,  carried  out  by  the  P.I.  and  his  colleague  Prof.  Edgardo 
Ficocelli  Varracchio,  were  published  in  J.  Phys.  B  16,  1097  (1983)  and 

released  as  Technical  Report  QFT-3  uf  ti.is  project.  During  his  visit  at 
the  university  of  Oregon  for  the  purpose  of  this  research  collaboration, 
Prof.  Ficocelli  Varracchio  also  wrote  a  review  on  the  f ield-theoretic 
approach  to  the  theory  of  three -body  atomic  collision  phenomena  based  on 
our  novel  second  quantization  representation.  This  was  published  in  Ann. 
Phys.  (N.Y.)  145 .  131  (1983)  and  released  as  Technical  Report  QFT-4  of  this 
contract.  The  same  general  approach  was  used  by  the  P.I.  and  his 
Postdoctoral  Research  Associate  Dr.  J.D.  Gilbert  to  derive  the  Fock-Tani 
representation  Hamiltonian  for  electron-hydrogen  scattering.  This  was 
released  as  Technical  Report  QFT-5. 


As  a  basis  for  generalization  of  our  calculations  on  atomic  collision 
phenomena  so  as  to  include  perturbing  effects  of  the  environment  (for 
example,  a  plasma)  in  which  the  colliding  atoms  and  ions  are  immersed,  the 
P.I.  and  his  Postdoctoral  Research  Associate  Dr.  Charles  F.  Hart  developed 
a  general  technique  of  "Liouwellian  Green's  functions"  for  calculation  of 
dissipative  energy-shift  and  line -broadening  effects  on  atoms  in  a 
perturbing  medium.  This  research  was  published  as  a  series  of  three 
papers:  Phys.  Rev.  A  28,  1056  (1983),  Phys.  Rev.  A  28,  1072  (1983),  and 
Phys.  Rev.  A  28,  186  (1983)  and  released  as  Technical  Reports  QFT-6,  QFT-7, 
and  QFT-8  of  this  contract.  An  abbreviated  account  of  this  work  was  also 
published  in  the  volume  Spectral  Line  Shapes .  Vol.  2  (Walter  de  Gruyter, 
Berlin,  1983)  which  constitutes  the  proceedings  of  an  International 
Conference  on  Spectral  Line  Shapes  at  which  the  P.I.  described  this  work. 

In  the  course  of  investigating  the  theory  of  the  charge -exchange 
reaction  D+  +  H  -*  D  +  H+,  the  P.I.  discovered  a  technique  for  elimination 
of  the  coordinates  of  any  one  of  a  system  of  N  interacting  particles  by 
performing  a  canonical  transformation  to  an  accelerated  coordinate  system 
centered  on  that  particle.  This  was  published  in  Phys.  Rev.  A  28,  3635 
(1983)  and  released  as  Technical  Report  QFT-9. 


Realizing  that  a  truly  fundamental  treatment  of  the  theory  of  atomic 
and  ionic  collisions  occurring  in  a  perturbing  medium  must  come  to  grips 
with  the  problem  of  how  to  define  and  calculate  the  perturbed  bound  wave 
functions  of  the  colliding  species,  the  P.I.  and  his  Postdoctoral  Research 
Associate  Dr.  Charles  F.  Hart  derived  a  new  variational  principle  for 
determination  of  such  perturbed,  decaying  (finite  lifetime)  states.  This 
was  published  in  Phys.  Rev.  Lett.  51,  1725  (1983)  and  released  as  Technical 
Report  QFT-10.  This  was  also  published  in  abbreviated  form  in  Quantum 
Electrodynamics  and  Quantum  Optics,  ed.  A.O.  Barut  (Proceedings,  NATO 
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Advanced  Study  Institute,  Boulder,  Colorado,  1983;  Plenum,  New  York,  1984), 
pp .  323ff . ,  at  which  conference  the  P.I.  was  an  invited  speaker. 

Building  on  our  previous  calculations  for  the  positronium  formation 
reaction  and  our  subsequent  work  on  similar  three  body  problems,  the  P.I., 
his  Postdoctoral  Research  Associates  Dr.  P.C.  Ojha  and  Dr.  J.D.  Gilbert, 
and  his  Ph.D.  dissertation  student  Jack  C.  Straton  formulated  an  improved 
first-order  theory  of  the  charge  transfer  reaction  H+  +  H(ls)  -*  H(ls)  +  H+ 
and  used  it  to  calculate  differential  and  total  cross  sections  for  this 
reaction.  The  results  were  found  to  be  much  more  accurate  than  the 
conventional  first  Born  approximation  and  to  give  very  good  agreement  with 
experiment  for  the  total  cross  section  at  energies  greater  than  10  keV  and 
for  the  differential  cross  section  at  energies  25,  60,  and  125  keV  in  an 
angular  range  of  **»  1  mrad  about  the  forward  direction.  This  work  was 
published  in  Phys .  Rev.  A  .33,  112  (1986)  and  released  as  Technical  Report 
QFT-11.  The  first-order  approach  was  subsequently  extended  by  the  P.I.  to 
a  formulation  of  the  distorted-wave  method  in  this  representation.  It  was 
published  in  Phys.  Rev.  A  33,  905  (1986)  and  released  as  Technical  Report 
QFT-12. 

The  P.I.  was  awarded  a  (unsolicited)  Humboldt  Prize  (U.S.  Senior 
Scientist  Award)  by  the  Alexander  von  Humboldt  Stiftung  for  the  period 
1984-85  and  spent  three  three  to  four-month  periods  working  in  the  Max 
Planck  Institut  fur  Strahlenchemie ,  Mulheim/Ruhr ,  West  Germany  during  the 
years  1984-86.  These  visits  were  funded  entirely  by  the  A.v.H.  Stiftung 
and  the  Max  Planck  Society.  Therefore  the  publications  resulting  from  the 
research  collaboration  between  the  P.I.  and  his  colleagues  there  do  not 
appear  as  ONR  Technical  Reports.  However,  a  few  of  these  publications  are 
closely  related  to  the  research  program  of  this  ONR  contract  and  are 
therefore  appended  to  the  Index  of  Publications  with  a  note  that  they  were 
not  funded  by  the  ONR. 

The  P.I.  's  Ph.D.  dissertation  student  Jack  C.  Straton  completed  his 
Ph.D.  dissertation,  entitled  "Charge  transfer  In  a  3  ■*  2-body,  reduced  mass 
Fock-Tani  representation:  first  order  results  and  an  introduction  to 
higher  order  effects",  based  on  research  funded  by  this  ONR  contract,  and 
received  his  Ph.D.  degree  in  June  1986.  He  is  currently  in  the  second  year 
of  an  appointment  as  a  National  Research  Council  -  NASA  Postdoctoral 
Research  Associate  at  the  NASA  Goddard  Space  Flight  Center  in  Greenbelt, 
MD,  where  he  works  in  the  atomic  theory  group  of  A.  Temkin,  R.J.  Drachman, 
and  A.K.  Bhatia.  Prior  to  taking  up  his  duties  under  that  appointment,  he 
completed  two  papers  based  on  his  research.  These  were  published  in  Phys. 
Rev.  A  35,  2729  (1987)  and  Phys.  Rev.  A  35,  3725  (1987)  and  released  as 
Technical  Reports  QFT-13  and  QFT-14. 

In  the  course  of  generalizing  our  treatment  of  reactive  collisions  of 
atoms  to  include  perturbing  effects  of  the  environment  (for  example,  a 
plasma)  we  realized  that  the  standard  Lippmann- Schwinger  theory  of 
collisions  does  not  apply  to  reactions  involving  metastable  or  unstable 
quasibound  states,  as  occur,  for  example,  in  molecular  predissociation, 
atomic  autoionization,  the  Auger  effect,  and  atomic  recombination  to  highly 
perturbed  and  broadened  atomic  states  in  a  plasma.  A  generalization  of 
Lippmann- Schwinger  collision  theory  applicable  to  such  cases  was  worked 
out,  published  in  Int.  J.  Quantum  Chem.  ££,  1483  (1986),  and  released  as 
Technical  Report  QFT-15. 
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The  calculations  on  atomic  charge -transfer  reactions  thus  far 
described  were  based  on  an  improved  first-order  theory  obtained  by 
canonically  transforming  the  Hamiltonian  and  then  applying  the  first  Born 
approximation  in  the  new  representation.  One  of  the  P.I.'s  Ph.D. 
dissertation  students,  Chaomei  Hu,  has  extended  this  approach  so  as  to 
include  an  important  class  of  second- Born  contributions  and  performed 
calculations  of  improved  differential  cross-sections  for  charge  exchange. 
Specifically,  she  included  both  initial  and  final -state  elastic  scattering 
contributions  in  a  post-prior  symmetrized  Fock-Tani  second  Born  calculation 
of  differential  cross  sections  for  the  charge  -  exchange  reaction  H++H-*H+H+. 
In  addition  she  calculated  post-prior  symmetrized  Fock-Tani  first  Born 
differential  cross  sections  for  the  reactions  H+  +  H(ls)  -*  H(nfim)  +  H+  and 
e+  +  H(ls)  -*  Ps(nfim)  +  H+  for  the  cases  n  -  1  and  n  —  2,  A  -  0,1.  Her 
Ph.D.  dissertation  will  be  completed  this  year  and  the  ONR  will  be  provided 
with  copies  of  the  dissertation  and  all  resultant  publications,  with 
acknowledgement  of  the  ONR  support. 

The  P.I.  and  another  of  his  Ph.D.  dissertation  students,  Fernando  A. 
Gutierrez,  began  a  program  of  calculation  of  cross  sections  for  those 
reactive  processes  of  atoms  in  plasmas  which  are  strongly  influenced  by  the 
plasma  environment.  Specifically,  we  developed  the  theory  of  a  new 
mechanism  for  atomic  recombination  to  highly-excited  and  perturbed  states 
of  hydrogen  in  a  plasma,  whereby  the  binding  energy  is  carried  away  by  a 
plasmon  (quantized  collective  plasma  oscillation)  rather  than  a  photon 
(radiative  recombination)  or  a  third  particle.  This  work  has  been 
submitted  to  Phys .  Rev.  A  and  a  preprint  is  appended  to  this  Report. 
Copies  of  the  dissertation  of  Mr.  Gutierrez  (expected  to  be  completed  this 
year)  and  any  resultant  publications  will  be  sent  to  the  ONR,  and  ONR 
support  of  this  work  will  be  acknowledged.  The  approach  generalizes  our 
earlier  work  on  reactive  scattering  by  large  systems  cited  as  item  17  in 
the  Index  of  Publications. 
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APPENDIX:  Work  on  plasmon-mediated  atomic  recombination 

in  progress  at  termination  of  the  contract 


Collective  mechanism  for  atomic  recombination  in  plasmas 


M.  D.  Girardeau  and  F.  A.  Gutierrez 

Department  of  Physics  and  Institutes  of  Chemical  Physics 

and  Theoretic?.!  Science 
University  of  Oregon,  Eugene,  Oregon  97403 

PACS  numbers:  34.90.+g,  34.60.+Z,  52.20.-j 
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Abstract 


A  new  mechanism  is  proposed  for  atomic  recombination  in 
plasmas,  whereby  the  binding  energy  is  carried  away  by  a  plasmon. 
It  is  suggested  that  this  mechanism  may  compete  with  radiative 
and  three-body  modes  in  the  case  of  recombination  to  sufficiently 
highly-excited  and  perturbed  states  of  hydrogen.  A  procedure  for 
calculation  of  the  transition  rate  is  outlined  in  a  model  which 
treats  the  plasma  oscillations  by  the  Bohm-Pines  canonical 
transformation  and  the  atomic  bound  states  by  a  second  canonical 
transformation . 


R 


Two  fundamental  quantities  in  a  classical  plasma  are  the 

2  % 

Debye  screening  length  XD=  (kgT/8rrne  )  and  the  (longitudinal) 
plasma  frequency1 


Uq  *  “p  (X  +  1  XDq2  )  '  KDq  <  1  (1> 

2  % 

where  <jp=(4nne  /ra)  .  The  first  term  in  (1)  arises  from  the  long- 
range  correlations  between  electrons  and  the  second  from  their 
thermal  motion-  We  consider  herein  only  the  simplest,  one- 
component  model;  accordingly,  e  and  m  are  the  electron  charge  and 
mass.  Quantization  gives  rise  to  plasmons,  the  longitudinal 
analog  of  the  usual  (transverse)  photons,  with  energies  - 

-11  y.  _3 

4x10  n  eV  where  n  ,  the  electron  number  density,  is  in  cm 


In  analogy  with  the  case  of  Cherenkov  photon  emission,  the 
condition  for  conservation  of  energy  and  momentum  in  the  emission 
of  a  real  plasmon  of  momentum  hq  and  energy  hw^  by  an  electron 
of  velocity  v  is2  u^=q-v  .  Since  there  is  a  maximum  plasmon  wave 
vector  '  '  t^iere  *-s  a  critical  velocity  vc'vup*D  below 

which  real  plasmon  emission  cannot  occur.  The  purpose  of  this 
note  is  to  suggest  and  briefly  analyze  a  new  collective  mechanism 
for  atomic  recombination  in  plasmas,  whereby  the  binding  energy 
is  carried  away  by  a  plasmon  rather  than  a  photon  or  a  third 
particle.  A  semiclassical  picture  of  the  process  is  shown  in  Fig. 
An  electron  in  an  unbound  (hyberbolic)  orbit  with  respect  to 
the  given  proton  approaches  it  with  a  small  enough  impact 
parameter  that  the  speed  of  the  electron  near  perihelion  exceeds 
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the  critical  velocity  vc  mentioned  above,  leading  to  emission  of 
a  real  plasmon  and  capture  of  the  electron  into  an  elliptic 
orbit. 

18  3 

Consider  a  hydrogen  plasma  with  n-'-lO  /cm  and  k^T~l  eV  ; 

then  \D~100  aQ  (aQ=Bohr  radius  of  the  unperturbed  hydrogen  atom) 
-2 

and  hOp^xlO  eV.  In  zeroth  approximation  one  may  determine  the 
perturbed  hydrogen  levels  and  wave  functions  by  solution  of  the 

Schrodinger  equation  with  the  screened  (Debye-Hiickel)  potential. 

4  5  .  —2 

One  finds  '  that  the  8s  level  has  binding  energy  -  4x10  eV  and 

the  9h  level  (the  highest  screened  Coulomb  bound  state  for  the 

-3  . 

given  vaiue  of  kD)  ^  8x10  eV  .  These  binding  energies  are  an 
order  of  magnitude  smaller  than  those  of  the  unperturbed  states, 
so  that  plasmon-mediated  recombination  to  these  levels  is 
energetically  allowed. 

A  correct  calculation  of  the  recombination  cross  section 
requires  a  quantum-mechanical  treatment  of  at  least  the 
electrons.  We  shall  use  a  simplified  model  consisting  of  a 
single,  fixed  proton  at  the  origin  immersed  in  a  medium  composed 
of  the  electrons  and  a  uniform,  immobile  positive  background 
which  crudely  models  the  other  protons  ("jellium") .  The 
interactions  of  the  electrons  with  each  other,  with  the  one 
chosen  proton,  and  with  the  rigid  positive  background  gives  rise 
to  collective  plasma  modes  which  are  crucial  for  this  mechanism 

A  h 

of  recombination.  The  Hamiltonian  H  in  Bohm-Pines  representation 
further  transformed  into  Fock  space  is 


in 


4  A  A  f 

hk  ,  c  and  c  annihilate  and  create  plasmons  with  momentum 

q  q 

hq  and  energy  hu  ,  q  is  a  wave  vector  cutoff  ~  X. \  n  is  the 
volume  of  the  system  (periodic  boundary  conditions) ,  N=nn  is  the 
number  of  electrons,  and 


fqkk'  =  (tJq  "  */m)  1  "  ["q  “  (V«)q*  &  +  &' )  ]~X  (3) 


In  (3)  we  have  neglected  terms  of  order  q  for  q<qc  ;  these 


represent  pure  quantum  effects  and  are  very  small  (less  than  1% 


of  <j  for  n  ~  10 

q 
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and  kT  ~  1  eV) .  In  (2)  represents  the  free 


electron  kinetic  energy,  Hpl  the  free  plasmons,  Ve_e  a  screened 

A 

electron-electron  interaction  and  Ve_p  the  corresponding  screened 


electron-proton  interaction,  Ve_pl  the  electron-plasmon 
interaction  representing  emission  and  absorption  of  single 


plasmons  by  electrons  in  the  field  of  the  proton,  V  _  -  the 

p  pi 

absorption  and  emission  of  virtual  plasmons  by  the  proton,  and  A 
a  constant  negative  energy  shift  which  results  from  the  long- 
range  correlations  between  the  electrons6. 


A  few  remarks  are  in  order: 

(a)  If  we  had  considered  a  proton  at  r  with  velocity  vq>v  , 

A  t  4  4 

then  Vp_pi  »  containing  a  factor  exp(-iq*rQ)  ,  would  describe 
Cherenkov-like  emission  of  real  plasmons  by  the  fast  proton.  (In 
fact,  this  was  the  original  motivation  of  Pines  ) .  In  the  present 
case  the  proton  is  at  rest  and  the  proton  emits  and  absorbs  only 
virtual  plasmons,  which  will  be  found  to  give  a  negative  energy 
shift. 
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3  A 

(b)  In  the  original  work  of  Pines  ve_p^  was  neglected  compared 
to  the  leading  term  giving  real  plasmon  emission  by  the  proton. 

A 

In  the  present  case,  since  vp_pi  gives  only  a  constant  energy 

A 

shift,  ve_pi  ^e  leading  term  for  the  dynamics  of  real  plasmon 
emission  by  electrons  in  the  field  of  the  proton. 

A 

(c)  In  coordinate  space  Ve_^  can  be  represented  to  a  good 

3  7 

approximation  '  by  the  potential 


(4) 


where  r^  is  the  distance  between  the  proton  and  the  ith  electron. 

(d)  The  derivation  of  (2)  and  (3)  is  rather  lengthy  and  is  given 
in  the  Bohm-Pines  papers  and  a  forthcoming  paper  on  the  details 
of  the  plasmonic  recombination  mechanism. 

A 

The  demonstration  that  vp_pl  gives  rise  to  a  negative  energy 
shift  proceeds  as  in  other  Bose  quasiparticle  theories,  by 
carrying  out  a  canonical  transformation  to  cancel  the  terms 

A  AAA 

Vp_pi  linear  in  plasmon  operators  cq  and  c  1  .  The  required 
(unitary)  transformation  is 


A  _  1  A  A  A  /  2  2  \SC  A_1  A  J.  A  A_1  A  A  * 

S  cqS  =  cq  *  (  2ireVnqZh«q  )  *  ,  S  1  c*S  =  (S  1  cqS)  f 


(5) 


and  leads  to  a  negative  shift  A  identical  with  (2): 


A_1  A 


A  A 


s~l<V  +  Vpi)s  “  +  1 


Pi 


(6) 
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As  in  similar  fixed-source  theories,  this  shift  is  interpreted  as 
the  energy  of  virtual  zero-point  plasmons  absorbed  and  emitted  by 
the  proton.  The  only  other  term  in  the  Hamiltonian  (2)  affected 


by  the  transformation  (5)  is  the  term  Ve_p^  which  is  transformed 

into  V  ,+V  '  where 

e-pl  e-p 


which  can  be  regarded  as  a  dynamical  correction  to  the  static 

screened  Coulomb  proton-electron  interaction  .  Inserting  the 

explicit  expression  (3)  ,  one  sees  that  the  summand  tends  to 

cancel  between  regions  with  opposite  alignments  of  the  relevant 

vectors,  whereas  (4ne2/q2)  in  V  is  isotropic.  It  is  therefore 

e~P 

within  the  spirit  of  the  RPA  to  neglect  V  _  '  compared  to  V  , 

e  p  e  p 

and  we  shall  do  so  herein.  Indeed,  if  one  neglects  q  compared  to 
k'  (since  q«Jc  #  k'  >qc)  and  neglects  (hq2/2m)  compared  to  cj^  as 
was  done  in  (3),  then  for  fixed  k  and  Jc'  the  summand  is  an  odd 

,  -4  A 

function  of  q  ,  so  that  ve_p'  vanishes  in  this  approximation. 

At  this  point  it  is  interesting  to  compare  our  result  with 
the  Ecker-Weizel  potential8'9 


VEW<r> 


-r/xD 

[  (l/r)e  u  +  (1/XD)J 


(7) 


(8) 


between  electrons  and  protons  In  a  plasma,  which  has  been  used  to 
interpret  the  plasma  shifts  of  the  discrete  and  continuous 
spectra  of  hydrogen.  If  we  use  the  Debye-Huckel  potential  for  the 

A 

screened  interaction  V  in  coordinate  space  and  add  the 

e-p 

negative  shift  2A  [one  A  from  the  proton  energy  shift.  Eg.  (6) , 
and  another  A  from  the  energy  shift  of  each  electron,  Eqs.  (1) 
and  (2) ]  then  we  obtain  an  effective  potential 

Veff(r)  =  _e  f<Vr)e  u  +  (0.7AD)]  (9) 

which  is  very  similar  to  (8). 

Another  change  of  representation  is  useful  for  extracting 

A 

from  Ve_pl  the  term  representing  free-bound  electron  transitions 
with  plasmon  emission.  The  required  canonical  transformation  in 
the  simplest  case  where  the  bound  state  is  a  single-particle 

state  (the  proton  being  treated  here  as  a  fixed  force  center)  is 

•  A  1011 
effected  by  a  unitary  operator  U  given  by  ' 


F 


A 


(10) 


A  A  | 

where  Ay  and  Ay  are  the  (Fermi)  annihilation  and  creation 
operators  for  electrons  bound  in  (here  perturbed)  hydrogen 
orbitals  centered  on  the  proton  at  the  origin: 


(11) 


K  - 1 


u  t 


The  index  v  stands  for  (nlm)  (and  the  spin  z-component  s,  which 
will  be  suppressed  herein) ;  $v  is  the  (perturbed)  hydrogen 
momentum  wave  function,  the  Fourier  transform  of  the  perturbed 
spatial  wave  function  4>v(x)  •  The  new  bound'  state  Fermi 

A  A  1 

annihilation  and  creation  operators  a^  and  ay  introduced  by  this 
transformation  anticommute  with  the  original  (plane-wave) 

A  AX  A  A  X 

electron  operators  e^  and  (although  the  e^  and  do  not 

A 

anticommute)  and  they  commute  with  the  plasmon  operators  and 

A  1 

Cg  .  The  canonical  transformation  effected  by  (10)  is 


A.*)  AAAr*  A  V-v  A 

u  eku  -  ek  -  2  kk,  V-  v  +  2  „^<k)  % 


(12) 


where  is  the  bound  state  kernel 


4k*-=2  K<k>  K<k'> 


(13) 


Insertion  of  (12)  into  Ve_pl  ,  Eq.  (2),  yields  several  terms;  the 
one  describing  the  plasmonic  recombination  and  the  inverse 
ionization  process  is 
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(14) 


V(e 


a  +  pi) 


-  I  (  Sq  K  (qv,V,k) 

(53«c) 


A 


+ 


h.c 


with 


(qv  |  V  |  k)  =  «  3/2  (  — T-5  )  1  (  ^2  )  [  K(k+k'"q)  *qkk' 


k'  >qr 


(k')' 


“5  ?l/(k"+k' ~q)  Cqk"k'  Ak'k  ]  (15) 


The  effect  of  the  subtraction  term  involving  A^,  ^  is  to 
orthogonal ize  the  k  dependence  of  the  matrix  element  to  the  bound 

A 

electron  subspace.  This  is  physically  correct  since  e^ 
annihilates  an  unbound  electron.  An  equivalent  point  of  view  is 
to  regard  e^  as  the  annihilation  operator  for  an  orthogonalized 
plane  wave  rather  than  a  (Born  approximation)  free  plane  wave. 


The  physics  of  the  emission  in  the  matrix  element  (15)  is 
contained  in  •  Interpreting  v^=ftj c/m  as  the  random 

(classical)  velocity  of  the  electrons  in  the  plasma,  one  sees 
that  the  first  term  in  (3)  has  a  singularity  at  •  This  is 

the  previously-mentioned  condition  for  real  plasmon  emission. 
Because  of  the  constraint  q<q  only  the  few  electrons  in  the 

C 
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high-velocity  tail  of  the  Maxwellian  distribution  will  be  able  to 
satisfy  this  condition.  For  this  reason  Bohm  and  Gross1  neglected 
the  singularity  in  getting  the  dispersion  relation  (1);  we 
continue  using  this  approximation  here.  Things  are  different  for 
the  second  term  in  (3),  where  the  singularity  occurs  at 

<Jq  =  q-  [  vk  +  (hic'/m)  ]  U6) 

The  term  hie'/®  represents  the  increase  in  velocity  of  the 
electrons  produced  by  the  short  range  (  k'  >qc  )  attractive 
interaction  with  the  proton;  it  is  crucial  for  the  mechanism  in 
which  we  are  interested,  because  now  electrons  with  any  velocity 
in  the  field  of  the  proton  are  potentially  able  to  produce  real 
plasmon  emission.  The  reason  is  that  the  sum  on  k'  appearing  in 
(15)  ensures  that  for  any  v^  there  is  always  a  k'  such  that  (16) 
is  satisfied.  This  is  consistent  with  the  semiclassical  picture 
of  Fig.  1,  according  to  which  an  electron  starting  with  a  low 
velocity  v^  can  increase  its  velocity  until  it  emits  a  plasmon 
and  becomes  bound  to  the  proton. 

In  zeroth  approximation  the  wave  functions  ^(k)  in  (15)  can 
be  taken  to  be  the  eigenstates  of  the  screened  Coulomb  potential 
satisfying 

(  h2k2/2m  )  ^(ic)  -  n"1  J  (  4ne2/q2  j  ^(ic+q)  -  ^(ic) 

q>qc 


(17) 


With  the  above  choice  one  neglects  the  imaginary  part  of  the 

atomic  self  energy,  i.e.,  the  finite  lifetime  of  the  perturbed 

atomic  state.  In  particular,  the  simplified  model  considered 

herein  neglects  the  Stark  broadening  due  to  the  other  positive 
12 

ions  since  they  are  treated  as  a  uniform  background.  The  Stark 

broadening  due  to  the  electrons  is  in  principle  included  if  one 

uses  the  appropriate  generalized  Schrodinger  equation  and  this 

broadening  is  expected  to  be  comparable  to  that  due  to  the 

positive  ions  at  the  densities  and  temperatures  considered  herein 
18  3 

(n  ~  10  /cm  ,  kfiT  ~  1  eV)  .  We  intend  to  generalize  the 
treatment  to  a  two-component  plasma  (both  protons  and  electrons 
treated  dynamically)  in  subsequent  work. 

The  leading  approximation  to  the  transition  rate  for 
plasmon-mediated  recombination  is 

W(qi>  |  k)  =  (2rr/h2)  |  (qviTlk)  |2<Nk(l  +  Nq)  (1  -  Ny)>  «(efi)  (18) 

where 

€fi  =  %  +  6„  “  (h2*2/2m) 


The  expression  (18)  is  the  product  of  the  standard  expression  for 
the  transition  rate  of  a  reaction  in  vacuum  by  a  statistical 
factor  <•••>  expressing  the  effects  of  occupancy  of  the  initial 
and  final  states.  In  thermal  equilibrium 
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(20) 


A 

<55 

V 

u 

A 

V 

1  +  <N  >  ) 

f  1  -  <N  >  ) 

*  \ 

q  / 

l  V  ) 

This  is  the  standard  statistical  factor  for  a  transition  in  a 
medium.  Since  the  model  only  allows  formation  of  a  single  atom 
(only  one  proton  at  the  origin,  all  others  a  "smeared 

A 

background")  and  <Nj/>  is  the  mean  occupation  of  level  v  before 

A 

formation  of  an  atom  in  this  state,  we  may  take  <Ny>=0  and 

A 

replace  the  Fermi  blocking  factor  (1  -  <N^>)  by  unity.  The  mean 

A 

number  <N^>  of  plasmons  at  each  given  q  is  also  «  1  in  the 
density-temperature  regime  under  consideration  (nondegenerate 
plasmon  gas) .  Then  (18)  reduces  to 

W (qy  | k)  %  (27r/h2)  |  (qi/|T|k)  |2fk  5(efi)  (21) 

,  A 

where  f^  is  the  electron  distribution  function  <Nk>  of  the 
plasma,  nearly  Maxwellian  in  the  given  regime.  In  leading  order 
the  T  matrix  (qviTik)  reduces  to  the  Hamiltonian  matrix  element 
(qi/|V|k)  for  the  given  reaction.  Although  this  has  the  appearance 
of  a  Born  approximation,  it  is  expected  to  be  much  better  than 
the  usual  "bare  Born"  approximation.  In  the  first  place,  the 
atomic  wave  function  is  the  screened  Coulomb  one,  which 
takes  into  account  in  zeroth  order  the  strong  perturbation  due  to 
the  plasma  environment.  In  the  second  place  the  matrix  element 
(15)  includes  the  orthogonalization  correction  which  is  neglected 
in  the  usual  Born  approximation.  A  similar  orthogonalization  term 
has  been  proposed  previously  as  a  correction  to  Born- 
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•  14 

approximation  photoionization  matrix  elements  .  In  the  present 
formulation  this  term  arises  naturally  as  a  consequence  of  the 
canonical  transformation  (12). 

The  numerical  evaluation  of  the  plasmonic  recombination  rate 
and  a  comparison  with  the  radiative  rate  is  underway;  this 
together  with  a  more  detailed  description  of  the  representations 
used  herein  will  be  given  in  a  subsequent  publication.  Ve  also 
plan  to  generalize  the  theory  to  the  more  realistic  two-component 
plasma  case. 
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Figure  Caption 


Fig.  l.  Semiclassical  picture  of  recombination  by  plasmon 

emission.  The  wavy  line  denotes  the  plasmon  and  the 
dashed  line  the  "post-collision"  hyperbolic  orbit  which 
would  have  been  followed  had  the  plasmon  not  been 
emitted. 
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